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Abstract. Star forming galaxies exhibit a variety of physical conditions, from 
quiescent normal spirals to the most powerful dusty starbursts. In order to study 
these complex systems, we need a suitable tool to analyze the information coming 
from observations at all wavelengths. We present a new spectro-photometric model 
which considers in a consistent way starlight as reprocessed by gas and dust. We 
discuss preliminary results to interpret some observed properties of VLIRGs. 

in 

o\ 

1. Introduction: why we care about nebular emission 

A common property of star forming galaxies is the presence of emission 
lines in the spectrum. Emission lines are produced by gas that is ion- 
Q H i ized by UV radiation emitted from young massive stars. Star- forming 

regions are the sites where the nebular emission is produced, but these 
regions are also dusty environments. It seems a natural goal for a model 
of star-forming galaxies to treat contextually the nebular emission and 
dust processing. 

We present the implementation of nebular emission calculation in 
the spectrophotometric code GRASIL (Silva et al., 1998). GRASIL is a 
code built to simulate the SED of galaxies including a careful treatment 
of starlight dust reprocessing. 

The study of emission lines with GRASIL will give us many further 
constraints about: 

— The dust obscuration: the optical thickness of star-forming regions 
can be derived from line ratio while it is very difficult to constrain 
optical thickness of ISM only from continuum. 

— The star formation rate: emission lines are believed to be powerful 
SFR estimators. We can now compare SFR derived from nebular 
lines with the results from FIR, UV etc. 
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2. The model 

In our model, the galaxy is composed by two main components: the 
bulge and the disk. Dust within the disk is divided in: i) molecular 
clouds (MCs), where the star formation is active, ii) diffuse medium 
(or cirrus). 

Young stars are supposed to be born into the MCs, and leave them 
progressively as their age increases. This is modelled by considering 
that a suitable fraction of the light of SSP is radiated inside the MCs. 
This fraction is a function of age, parametrized by an escaping time. 

So we have two populations: A) young stars inside MCs, B) older 
stars outside MCs. 

A molecular cloud is modelled as a thick spherical shell of dense 
gas (and dust) around a central point source representing all the star 
content of the cloud. The MCs, as well as old stars, are embedded in 
the diffuse medium. For more details on energy transport, refer to Silva 
et al. 1998 or Silva et al., this conference. 

To calculate emission lines we used the photoionization code CLOUDY 
94 (Ferland, 1996). 

We built a library of HII regions in order to avoid the computation of 
entire nebular emission via CLOUDY at every run of the code GRASIL. 

The HII regions of the library are calculated for different density, 
metallicity and radii. The SEDs of ionizing sources of HII regions are 
obtained by integrating the light of each of the two populations for 
several assumption of escaping time and IMF. 

The SEDs of ionizing sources are parametrized by the number of ion- 
izing photons for HI, Hel and Oil (Q(H), Q(He), and Q(O)). Q(H) gives 
the mass of ionized gas, while the ratios Q(He)/Q(H) and Q(0)/Q(H) 
are related to the hardness of ionizing flux and to the degree of ioniza- 
tion. In table I it is reported the list of computed lines. 

GRASIL calculates the Qs of each population (inside and outside 
the MC), interpolates on the library of HII regions, and computes the 
nebular emission for each population. Then, the nebular emission is 
extincted in the same way as the population that produced it. 

3. An application to Very Luminous InfraRed Galaxies 

VLIRGs (Lir > 1O 11,5 L0) are the most powerful star forming galaxies 
at low redshift. Recent results show that: VLIRGs do not follow the 
Meurer's relation Lfir/Xfuv versus UV spectral index (3 (Meurer et 
al., 1999; Meurer et al., 2000); SFR computed from Ha, even corrected 
for extinction, is always lower than SFR derived from FIR (Poggianti 
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Table I. Computed lines: Hydrogen recombination lines (upper pannel), 
Elium and metal lines (lower pannel) 



Lyct 1216 


Ly/3 1025 


Ly 7 972 


Ly5 949 


Ly 937 


Ly 930 


Ly 926 


Ly 922 


Ha6563 


H/3 4861 


H 7 4340 


HiS 4102 


H 3970 


H 3889 


H 3835 


H 3798 


Paci 18752 


Pa,3 12819 


Pa 7 10939 


Pa.5 10050 


Pa 9546 


Pa 9229 


Pa 9015 


Pa 8863 


Bret 40515 


Br/3 26254 


Br 7 21657 


Bri 19447 


Br 18175 


Br 17363 


Br 16808 


Br 16408 


Pfa 74585 


Pf/3 46529 


Pf 7 37398 


Pf5 32964 


Pf 30386 


Pf 28724 


Pf 27577 


Pf 26746 


Hio 123690 


Hu/3 75011 


n u 7 Ow* i 


Pin X ^1 977 


Tin 4fi71fi 


xx u / t)U 


u„ /ii 7nn 


xxu /±u.iU± 


Hell 1640 


Hell 1217 


HoII 1085 


Hell 4686 


Hell 3203 


Hell 2733 


Hell 2511 


Hcl 4471 


Hoi 5876 


Hcl 6678 


Hcl 10830 


Hel 3889 


Hcl 7065 


[CIJ9850 


[CI] 8727 


[CI] 4621 


[CI]609pm 


[CI]369pm 


[CII] 157.7pm 


CII] 2326 


CIII] 1908 


[NI]5199 


[NI]3466 


[NI] 10400 


[Nil] 6584 


[Nil] 6548 


[NII]5755 


[Nil] 122pm 


[NII]205pm 


NII]2141 


[NIII]57pm 


[OI]6300 


[OI]6363 


[OI]5577 


[OI]63pm 


[OI] 145pm 


[OIIJ3727 


[OH] 7325 


[OII]2471 


OIII]1663 


[OIII]5007 


[OIII]4959 


[OIII]4363 


[OIII]2321 


[OIII]88pm 


[0111152pm 


[OIV]26pm 


[Nell] 13pm 


[Nelll] 15.5pm 


[NeIII]36pm 


[NeIII]3869 


[NcIII]3967 


[NeIII]3343 


[NeIII]1815 


[NeIV]2424 


[NeIV]4720 


MgII2800 


[SiII]35pm 


[SII] 10330 


[SII]6731 


[SII]6717 


[SII]4070 


[SIIJ4078 


[SIII]19 


[SHI] 33. 5 


[SIII]9532 


[SIII]9069 


[SIII]6312 


[SHI] 3722 


[SIV] 10.4pm 


[ArII]69850 


[ArIII]7135 


[ArIII]7751 


[ArIII]5192 


[ArIII]3109 


[ArIII]3005 


[ArIII]22pm 


[ArIII]9pm 



& Wu, 2000). Can the above problems be solved with a correct picture 
of the obscuration? 

We have simulated with our code a starburst galaxy (barionic mass 
of 5 • 1O 1O M0) with a quiescent star formation lasting all its history, 
plus a final (analitycal) burst. We supposed that, as the burst begins, 
half of total (1.2 • 1O 1O M ) gas is under molecular form, and a fraction 
(from 70% to 99%) of it is converted in stars during the burst. As 
the molecular gas is consumed during the burst, the MCs become more 
and more transparent. This is to mimic the consumption of the gas and 
the feedback of SNe. The initial optical thickness of MC is chosen to 
reproduce observed eq. widht of Ha+[NII] and H/3 (-62.5 A and 0.7 A) 
and the ratio Lfir/^v (88, Poggianti et al., 2001). The time evolution 
of the model is shown in figure 1. 

The temporal evolution in the /3-FIR/FUV plane can be described 
in three phases: 

1) VLIRG phase: the burst is very obscured; IR emission comes from the 
burst while UV comes from the disk. There is no correlation between 
UV and IR because they come from different component of the galaxy. 

2) UV-bright phase: the burst has almost consumed the gas and the 
MCs become transparent; the model moves towards Meurer's relation. 

3) Final phase: the burst is aging and the model departs from the 
relation. The galaxy can not be observationally selected as a starburst 
any longer. 
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Figure 1. Left: The time evolution of the model in /3-FIR/FUV plane. Observed data 
are ULIRG observations (filed squares, Meurer et al. 2000) and UV-selected starburst 
(open circles, Meurer et al. 2000). Different paths correspond to different fraction 
of gas converted in stars. Right: SFR derived from Ha corrected for extinction with 
Balmer decrement versus SFR from FIR. 

In inital VLIRG phase, Ha corrected for extinction from Balmer 
decrement Ha/H/3 can be a wrong estimate of the "true" Ha (and then 
of the SFR, see fig. 1 right). In fact in this phase H/3 can be dominated 
by the population outside MCs affected by a low extinction; then, the 
Balmer decrement is altered because of the selective extinction. This 
effect can explain the low Ln ao /LBoi observed in VLIRGs (Meurer & 
Seibert, 2001). Concluding, we are able to interpret different UV, opti- 
cal, IR properties of obscured galaxies within a unique star formation 
selective (in age) extinction scenario. We found that the relation /3 - 
FIR/FUV is not unique in the obscured phase. For the same reason, the 
correction of Ha from balmer decrement can be wrong in this phase. 
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